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Objective:  To  improve  our  ability  to  identify  physiologic  dete¬ 
rioration  caused  by  critical  illness,  we  applied  nonlinear  and 
frequency-domain  analytical  methods  to  R-to-R  interval  (RRI)  and 
systolic  arterial  pressure  (SAP)  time  series  during  hemorrhagic 
shock. 

Design:  Prospective,  randomized,  controlled  trial. 

Setting:  Animal  laboratory  of  a  government  research  institute. 

Subjects:  Twenty  swine  (weight  36.4  ±  0.11  kg). 

Interventions:  Fixed-volume  hemorrhage  followed  by  resusci¬ 
tation;  off-line  analysis  of  RRI  and  SAP  data. 

Measurements  and  Main  Results:  Anesthetized  swine  (shock 
group,  n  =  12)  underwent  withdrawal  of  30  mL/kg  blood  in  10 
mL/kg  decrements.  A  control  group  (n  =  8)  received  maintenance 
fluids  only.  Electrocardiogram  and  arterial  pressure  waveforms 
were  acquired  at  500  Hz.  Eight  hundred-beat  data  sets  were 
analyzed  at  six  time  points:  at  baseline,  after  each  blood  with¬ 
drawal,  after  lactated  Ringer’s  resuscitation,  and  after  infusion  of 
shed  blood.  Nonlinear  methods  were  used  to  estimate  the  com¬ 
plexity  (approximate  entropy,  sample  entropy,  Lempel-Ziv  en¬ 
tropy,  normalized  entropy  of  symbol  dynamics),  RRI  bits  per  word, 


and  fractal  dimension  by  curve  lengths  and  by  dispersion  analysis 
of  the  RRI  and  SAP  time  series.  Fast  Fourier  transformation  was 
used  to  measure  the  high-frequency  and  low-frequency  powers  of 
RRI  and  SAP.  Baroreflex  sensitivity  was  assessed  in  the  time 
domain  with  the  sequence  method.  Hemorrhagic  shock  caused 
decreases  in  RRI  complexity  as  quantified  by  approximate  en¬ 
tropy,  sample  entropy,  and  symbol  dynamics;  these  changes  were 
reversed  by  resuscitation.  Similar  but  statistically  insignificant 
changes  in  fractal  dimension  by  curve  lengths  were  seen.  RRI 
high-frequency  power  decreased  with  hemorrhagic  shock — indi¬ 
cating  withdrawal  of  vagal  cardiac  input — and  was  restored  by 
resuscitation.  Similar  changes  in  baroreflex  sensitivity  were  seen. 
Hemorrhagic  shock  did  not  affect  SAP  complexity. 

Conclusions:  Hemorrhagic  shock  caused  a  reversible  decrease 
in  RRI  complexity;  these  changes  may  be  mediated  by  changes  in 
vagal  cardiac  control.  Assessment  of  RRI  complexity  may  permit 
identification  of  casualties  with  hemorrhagic  shock.  (Crit  Care 
Med  2007;  35:519-525) 

Key  Words:  shock;  hemorrhage;  electrocardiography;  heart  rate 
variability;  fractals;  nonlinear  dynamics;  spectrum  analysis 


Noninvasive  monitoring  may 
revolutionize  the  care  of  the 
critically  ill  if  data  analysis 
techniques  can  be  validated 
that  identify  physiologic  deterioration  be¬ 
fore  changes  in  conventional  vital  signs. 


‘See  also  p.  656. 
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In  pursuit  of  this  goal,  frequency-domain 
measures  of  heart  rate  variability  (HRV) 
have  been  proposed,  enabling  noninva¬ 
sive  quantification  of  the  autonomic 
modulation  of  the  heart  rate  (1,  2).  HRV 
metrics  have  been  associated  with 
changes  in  volume  status  (3, 4)  and  in  the 
depth  of  anesthesia  (5),  “inappropriate” 
vagal  predominance  has  been  associated 
with  adverse  outcomes  in  trauma  pa¬ 
tients  (6),  and  loss  of  HRV  has  been  pre¬ 
dictive  of  death  in  critically  ill  humans  (7, 
8).  The  frequency-domain  approach  has 
also  been  criticized  (9).  The  results  ob¬ 
tained  are  affected  by  the  methodology 
used  (10);  data  analysis  assumes  signal 
stationarity,  a  condition  difficult  to 
achieve;  and  the  approach  has  not  been 
designed  to  unveil  the  inherent  dy¬ 
namic  properties  of  complex  cardiovas¬ 
cular  regulatory  systems  (11).  As 
pointed  out  by  Goldberger  and  West 
(12)  and  Buchman  (13),  physiologic 
systems  exhibit  marked  signal  variabil¬ 
ity,  constantly  change  over  time,  and 
operate  in  nonlinear  states. 


In  linear  systems,  the  magnitude  of 
response  is  proportionate  to  the  strength 
of  the  stimulus.  In  a  nonlinear  system,  by 
contrast,  proportionality  does  not  hold; 
the  components  of  this  system  are  cou¬ 
pled  and  interact,  causing  abrupt  changes 
(14).  Such  systems  are  better  described 
by  nonlinear  statistics  that  determine  the 
degree  of  regularity  or  irregularity  in  the 
system  (15). 

Because  frequency-domain  methods 
assume  system  linearity,  nonlinear  statis¬ 
tical  methods  have  increasingly  been  ap¬ 
plied  to  HRV  analysis  (5,  7,  11,  16-18). 
Use  of  these  methods  is  predicated  on  the 
concept  that  complexity  is  a  fundamental 
characteristic  of  normal  physiologic  data 
such  as  the  R-to-R  interval  (RRI)  time 
series  (16),  whereas  loss  thereof,  or  “de- 
complexification,”  is  a  feature  of  disease 
(12)  and  of  impaired  adaptation  to  phys¬ 
iologic  stress  (19).  For  example,  loss  of 
RRI  complexity  as  manifested  by  de¬ 
creases  in  approximate  entropy  (ApEn), 
and  alterations  in  fractal  scaling  proper¬ 
ties,  have  both  been  identified  with  aging 


Crit  Care  Med  2007  Vol.  35,  No.  2  519 

Copyright  ©  Lippincott  Williams  &  Wilkins.  Unauthorized  reproduction  of  this  article  is  prohibited 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

01  FEB  2007  N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Loss  of  complexity  characterizes  the  heart  rate  response  to  experimental 
hemorrhagic  shock  in  swine 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Batchinsky  A.  I.,  Cooke  W.  H.,  Kuusela  T.,  Cancio  L.  C., 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

United  States  Army  Institute  of  Surgical  Research,  JBSA  Fort  Sam 
Houston,  TX  78234 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  |J|J 

unclassified  unclassified  unclassified 

7 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


(19-21).  Decrease  of  RRI  fractal  scaling 
was  documented  in  heart  failure  and  in 
Chagas’  disease  (22);  was  predictive  of 
mortality  after  myocardial  infarction  (18, 
23)  and  stroke  (24),  as  well  as  of  ischemia 
after  coronary  artery  bypass  grafting  (25); 
and  distinguished  subjects  with  coronary 
artery  disease  from  healthy  controls  (17). 

Reports  on  the  assessment  of  com¬ 
plexity  during  hypovolemia  or  hemor¬ 
rhagic  shock  (HS)  are  scarce.  Decreased 
ApEn  following  prolonged  bed  rest  was 
associated  with  reduced  tolerance  to  cen¬ 
tral  hypovolemia  induced  by  lower  body 
negative  pressure  (LBNP)  (26).  Decrease 
in  the  complexity  of  the  HR  but  not  of  the 
blood  pressure  has  been  documented  in 
healthy  humans  during  orthostatic  stress 
secondary  to  LBNP  (3).  Hemodialysis- 
induced  hypotension  decreased  ApEn  in 
patients  with  and  without  diabetes  melli- 
tus  (27).  In  dogs,  Palazzolo  et  al.  (28) 
identified  a  decrease  in  RRI  ApEn  and 
symbol  dynamics  entropy  (SymDyn)  sec¬ 
ondary  to  hypotension. 

The  objective  of  this  study  was  to  de¬ 
termine  whether  loss  of  RRI  and  systolic 
arterial  pressure  (SAP)  complexity  is  a 
characteristic  of  HS.  We  hypothesized 
that  hemorrhage  causes  loss  of  complex¬ 
ity  in  these  signals  and  that  resuscitation 
restores  it.  Several  different,  but  poten¬ 
tially  complementary,  nonlinear  methods 
were  used. 

MATERIALS  AND  METHODS 

This  study  was  approved  by  the  U.S.  Army 
Institute  of  Surgical  Research  Animal  Care 
and  Use  Committee  and  was  carried  out  in 
accordance  with  the  guidelines  set  forth  by  the 
Animal  Welfare  Act,  other  federal  statutes  and 
regulations,  and  the  1996  Guide  for  the  Care 
and  Use  of  Laboratory  Animals  of  the  National 
Research  Council. 

Animal  Preparation.  Twenty  female  York¬ 
shire  pigs  (hemorrhage  group  n  =  12,  control 
group  n  =  8)  weighing  36.4  ±  0.11  kg  sem 
were  fasted  for  24  hrs,  anesthetized  with  a 
combination  of  isoflurane  (0.62  ±  0.04%)  and 
intravenous  ketamine  (332.1  ±  10.4  sem  pig / 
kg/hr)  (Ketamine  Hydrochloride,  St.  Joseph, 
MO),  and  volume-control  ventilated  with 
room  air  at  a  tidal  volume  of  10  mL/kg  and 
respiratory  rate  of  12  per  minute.  Respiratory 
rate  was  adjusted  at  baseline  to  maintain  nor¬ 
mocapnia  (Paco2  =  35-45  mm  Hg)  and  was 
kept  constant  throughout  the  study  to  stan¬ 
dardize  the  effect  of  respiration  on  HRV  (29, 
30).  Bilateral  surgical  catheterizations  of  the 
carotid  arteries  and  external  jugular  veins  and 
a  midline  laparotomy  and  splenectomy  (to 
preclude  autoresuscitation  during  hemor¬ 
rhage)  were  performed.  Each  animal  received 


a  bolus  of  lactated  Ringers’  solution  (LR) 
equivalent  to  1.5  times  the  spleen  weight.  Ar¬ 
terial  blood  pressure  (ABP)  was  measured  via 
the  carotid  line  and  was  transduced  (Transpac 
IV,  Abbott,  North  Chicago,  IL)  at  heart  level. 
The  heart  rate  was  obtained  from  the  electro¬ 
cardiogram. 

Experimental  Protocol.  After  1  hr  of  base¬ 
line  stabilization,  the  experimental  protocol 
was  started  (Table  1).  In  the  hemorrhage 
group,  blood  was  manually  withdrawn  via  a 
syringe  from  the  left  carotid  artery  catheter 
during  each  of  three  10-min  bleeding  periods 
at  a  constant  rate  of  1  mL/kg/min,  for  a  total 
bleed  of  30  mL/kg.  This  approximated  a  45% 
hemorrhage.  The  blood  was  collected  into 
commercial  collection  bags  containing  CPD 
anticoagulant  (Terumo,  Tokyo,  Japan).  Each 
bleed  period  was  followed  by  a  15-min  stabili¬ 
zation  period.  After  the  third  bleed,  the  pigs 
were  resuscitated  with  two  times  the  shed 
blood  volume  of  LR,  followed  by  reinfusion  of 
shed  blood.  LR  resuscitation  and  blood  rein¬ 
fusion  were  performed  with  the  245  Ranger 
(Blood/Fluid)  Warming  System  (Augustine 
Medical,  Eden  Prairie,  MN).  The  animals  were 
observed  for  30  mins  thereafter  and  then  were 
killed  by  an  overdose  of  sodium  pentobarbital 
(Fatal-Plus,  Dearborn,  MI).  Animals  in  the 
control  group  received  a  maintenance  LR  in¬ 
fusion  at  a  2  mL/kg/hr  rate  only. 

Waveform  Analysis.  Electrocardiogram  and 
ABP  waveforms  were  continuously  monitored 
throughout  the  experiment  and  recorded  at  500 
Hz  for  off-line  analysis  during  discrete  15-min 
observation  periods.  This  allowed  us  to  do  the 
analysis  on  relatively  stable  data  sets  collected 
during  periods  between  interventions  (i.e.,  epi¬ 
sodes  of  blood  withdrawal,  resuscitation,  or 
blood  reinfusion).  Six  time  points  (baseline,  after 
bleed  1,  after  bleed  2,  after  bleed  3,  after  LR 
resuscitation,  and  after  blood  reinfusion,  Table 
1)  were  analyzed  for  every  subject.  Next,  800- 
beat  data  sets  from  within  each  15-min  wave¬ 
form  recording  were  imported  into  WinCPRS 
software  (Absolute  Aliens  Oy,  Turku,  Finland). 
Eight  hundred  beats  were  selected  because  non¬ 
linear  statistics  are  sensitive  to  changes  in  the 


number  of  data  points  (especially  ApEn).  Auto¬ 
matic  identification  of  R-waves  was  carried  out 
by  means  of  an  isoelectric  line-shift  algorithm. 
We  manually  verified  correct  identification  of 
R-waves  by  the  software  in  every  data  set.  When¬ 
ever  possible,  ectopic  beats  (EBs)  were  avoided 
by  selecting  an  800-beat  section  free  from  ec- 
topy.  When  present,  EBs  were  replaced  by  a 
linear  interpolation  function  of  the  software. 
Of  the  20  animals  included  in  this  study, 
EBs  were  identified  at  a  few  time  points  in 
ten  of  them.  However,  in  most  cases  (eight 
of  ten),  only  one  EB  was  identified  within 
the  800-beat  segment  analyzed,  accounting 
for  a  0.13%  interpolation  rate  per  800  beats. 
In  the  remaining  two  of  ten  animals  with 
EBs,  three,  12,  and  15  EBs  were  identified 
within  three  separate  800-beat  segments,  ac¬ 
counting  for  0.38%,  1.5%,  and  1.8%  inter¬ 
polation  rates,  respectively. 

The  software  generated  the  instantaneous 
RRI  and  SAP  time  series.  Before  spectral  calcu¬ 
lations,  signals  were  windowed  with  Hanning 
function  and  resampled  with  a  five  sample  per 
second  rate.  The  resampled  signals  were  zero- 
padded  to  the  next  power  of  2  before  the  FFT 
routine  was  used  to  calculate  the  power  spectra. 
Next,  the  total  power  spectral  density  (TP,  fre¬ 
quency  range:  0.003-0.4  Hz),  the  low-frequency 
power  (LFP;  frequency  range,  0.04-0.15  Hz) 
and  high-frequency  power  (HFP;  frequency 
range,  0.15-0.4  Hz)  components  of  the  RRI  and 
SAP  spectra  were  calculated  for  each  data  set 
with  baseline  normalization.  The  LFP/HFP  ratio 
of  the  RRI  was  taken  as  an  index  of  sympatho- 
vagal  balance  (1,  2).  The  integrated  baroreflex 
sensitivity  (BRS),  ApEn,  sample  entropy 
(SampEn),  Lempel-Ziv  entropy,  SymDyn,  RRI 
bits  per  word,  fractal  dimension  by  curve  lengths 
(FDCL),  and  fractal  dimension  by  dispersion 
analysis  were  calculated  for  the  RRI  and  SAP 
time  series  in  the  same  800-point  data  sets  using 
the  WinCPRS  software,  as  described  by  Kuusela 
et  al.  (11)  (Appendix).  Before  entropy  calcula¬ 
tions,  linear  trends  were  removed. 

Statistical  Analysis.  Nonnormal  data  were 
logjQ-normalized  if  possible.  Independent 
samples  /-tests  or  Mann-Whitney  U  tests 


Table  1.  Experimental  protocol 


Duration, 

Time  Point  Mins  Event 


Baseline 

30 

Hemorrhage  1 

10 

Recovery  1 

15 

Hemorrhage  2 

10 

Recovery  2 

15 

Hemorrhage  3 

10 

Recovery  3 

15 

LR  infusion 

25 

Observation 

30 

Blood  reinfusion 

25 

Observation 

30 

Total  study  time 

215 

EKG  and  ABP  waveform  recording 
Constant  rate  1  ml-kg-min  bleed 
Observation/EKG  and  ABP  waveform  recording 
Constant  rate  1  ml-kg-min  bleed 
Observation/EKG  and  ABP  waveform  recording 
Constant  rate  1  ml-kg-min  bleed 
Observation/EKG  and  ABP  waveform  recording 
1.5  times  the  shed  volume  of  LR  infused 
Observation/EKG  and  ABP  waveform  recording 
Reinfusion  of  shed  blood 
Observation/EKG  and  ABP  waveform  recording 


EKG,  electrocardiogram;  ABP,  arterial  blood  pressure;  LR,  lactated  Ringers’  solution. 
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Table  2.  Study  results 


Group 

Baseline 

10  mL 

20  mL 

30  mL 

Post-LR 

Post-Tx 

p* 

jj' i*  -fc 

RRI  mean 

C 

618 

629 

634 

646 

649 

640 

B 

SAP  mean 

676 

600 

475* 

388** 

493*** 

570 

.002 

<.001 

.003 

C 

94 

97 

98 

98 

89 

88 

B 

RRI  ApEn 

100 

77* 

70** 

53** 

80 

105*** 

.002 

<.001 

.012 

C 

1.26 

1.20 

1.12 

1.04 

1.14 

1.08 

B 

RRI  SampEn 

1.21 

1.03 

1.00 

0.76* 

0.99 

1.09 

.017 

C 

1.35 

1.28 

1.22 

1.10 

1.20 

1.12 

B 

RRI  SymDyn 

1.29 

1.04 

0.99 

0.71* 

0.95** 

1.11 

.010 

.017 

C 

0.73 

0.72 

0.66 

0.67 

0.67 

0.67 

B 

RRI  FDCL 

0.68 

0.58* 

0.61 

0.50** 

0.64 

0.70 

.034 

.013 

C 

1.88 

1.86 

1.81 

1.77 

1.79 

1.76 

B 

RRI  BPW 

1.82 

1.74 

1.74 

1.62 

1.76 

1.79 

NS 

C 

4.40 

4.32 

3.99 

4.02 

4.04 

4.00 

B 

RRI  FDDA 

4.10 

3.49* 

3.69 

3.00** 

3.83 

4.20 

.034 

.013 

C 

1.15 

1.21 

1.20 

1.16 

1.12 

1.18 

B 

RRI  HFP 

1.07 

1.06* 

1.06** 

1.05*** 

1.05t 

1.08 

.012 

.025 

.008 

C 

8.36 

9.95 

9.58 

11.55 

9.69 

7.94 

B 

SAP  HFP 

3.76 

2.27 

1.02* 

0.23** 

2.21*** 

3.65 

.004 

<.001 

.015 

C 

4.50 

5.11 

5.03 

4.38 

5.47 

7.85 

B 

5.69 

9.90 

23.90* 

27.25** 

18.49*** 

2.84f 

.001 

.019 

.012 

BRS 

C 

1.33 

1.36 

1.39 

1.35 

1.26 

0.83 

B 

0.94 

0.61* 

0.24** 

0.15** 

0.36*** 

1.02 

.029 

<.001 

<.001 

Baseline,  baseline  data  collection;  10  mL,  20  mL,  30  mL,  data  collection  after  10  mL,  20  mL,  or  30  mL,  respectively,  of  blood  removed;  Post-LR,  after 
resuscitation  with  Ringers’  lactate;  Post-TX,  after  blood  reinfusion;  p*,  p**,  p***,  and  pt,  levels  of  significance  for  values  in  the  bleed  vs.  control  groups 
at  each  time  point  by  independent  samples  or  Mann  Whitney  U  tests;  RRI.  R-to-R  interval;  C,  control  group;  B,  bleed  group;  SAP,  systolic  arterial  pressure; 
ApEn,  approximate  entropy;  SampEn,  sample  entropy;  SymDyn,  symbol  dynamics  entropy;  FDCL,  fractal  dimension  by  curve  lengths;  BPW,  bit  per  word; 
FDDA,  fractal  dimension  by  dispersion  analysis;  HFP,  high-frequency  power;  BRS,  baroreflex  sensitivity. 

Data  are  mean  ±  sem. 


1  2  3  4  5  6 

Timepoint 

Figure  1.  Changes  in  R-to-R  interval  over  time.  Triangles,  hemorrhaged  animals;  circles,  controls.  For 
statistical  significance,  see  Table  2.  f-tests:  *p  =  .002;  **p  <  .001;  ***p  =  .003. 


were  performed  at  each  time  point,  control 
vs.  bleed  groups.  Significance  was  accepted 
at  p  <  .05. 

RESULTS 

Study  results  are  summarized  in  Ta¬ 
ble  2.  In  response  to  hemorrhage,  the 
RRI  decreased  linearly,  differed  from 
controls  after  the  second  bleed,  and  was 
restored  to  control  values  by  reinfusion 
of  shed  blood  (Table  2  and  Fig.  1).  SAP 
(Fig.  2)  differed  from  controls  after  the 
first  bleed  and  was  restored  by  LR  re¬ 
suscitation. 

The  ApEn,  SampEn,  and  SymDyn  (Fig. 
34  and  3B;  Table  2)  of  the  RRI  all  showed 
similar  graphic  appearance,  decreased  with 
hemorrhage,  and  were  restored  by  resusci¬ 
tation.  An  analogous  trend  was  observed  in 
FDCL  (Fig.  3 C),  although  these  differences 
were  not  significant.  RRI  bits  per  word  dif¬ 
fered  after  the  first  and  third  bleeds  (Table 


2).  The  values  of  fractal  dimension  by  dis¬ 
persion  analysis  were  systematically  lower 
in  the  hemorrhaged  animals  and  differed 
from  controls  after  the  first  bleed  (Table  2). 
RRI  Lempel-Ziv  entropy  did  not  change 
(data  not  shown).  Changes  in  SAP  ApEn, 
SampEn,  SymDyn,  Lempel-Ziv  entropy, 
FDCL,  fractal  dimension  by  dispersion 


analysis,  and  bits  per  word  were  not  signif¬ 
icant  and  are  not  shown. 

Spectral  analysis  revealed  a  significant 
decrease  in  RRI  HFP  (Fig.  44),  which  was 
reversed  by  resuscitation.  SAP  HFP  in¬ 
creased  with  hemorrhage,  differed  be¬ 
tween  groups  after  bleed  2,  and  re¬ 
mained  significantly  increased  until 
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Figure  2.  Changes  in  systolic  arterial  pressure  over  time.  Triangles,  hemorrhaged  animals;  circles, 
controls,  f-tests:  *p  =  .002;  **p  <  .001;  ***p  =  .012. 
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Figure  3.  A,  changes  in  R-to-R  interval  approximate  entropy  over  time.  Triangles,  hemorrhaged 
animals;  circles,  controls.  Mann-Whitney  U  tests:  *p  =  .017.  B,  changes  in  RRI  sample  entropy  over 
time.  Triangles,  hemorrhaged  animals;  circles,  controls.  Mann-Whitney  U  tests:  *p  =  .010;  **p  =  .017. 
C,  changes  in  RRI  fractal  dimention  by  curve  lengths  over  time.  Triangles,  hemorrhaged  animals; 
circles,  controls.  Not  significant. 


blood  reinfusion  (Fig.  45).  SAP  TP  in¬ 
creased  in  concert  with  the  SAP  HFP 
(data  not  shown).  Changes  in  RRI  LFP/ 
HFP  ratio,  RRI  LFP,  RRI  TP,  and  SAP 
LFP  were  not  significant  (data  not 
shown).  BRS  decreased  with  hemor¬ 
rhage  and  was  restored  by  resuscitation 
(Fig.  4C). 


DISCUSSION 

The  principal  finding  of  this  study  is 
that  RRI  complexity,  as  assessed  by  sev¬ 
eral  complementary  measures,  decreases 
after  controlled  30  mL/kg  hemorrhage  in 
anesthetized  and  mechanically  ventilated 
swine.  Loss  of  complexity  occurred  in  as¬ 


sociation  with  decreases  in  vagal  cardiac 
oscillations  (as  quantified  by  the  RRI 
HFP)  and  in  cardiovagal  baroreflex  sensi¬ 
tivity.  All  of  these  values  were  restored  by 
resuscitation. 

One  of  the  most  frequently  used  com¬ 
plexity  measures  was  introduced  by  Pin- 
cus  (31)  and  involves  quantification  of 
the  degree  of  entropy  in  the  system  as  an 
approximate  measure  of  complexity.  Hor- 
nero  et  al.  (32)  stated  that  ApEn  is  a 
robust  approximate  complexity  metric  of 
a  signal  that  combines  information  on 
HRV  and  the  number  of  signal  harmonics 
and  correlates  with  the  signal-noise 
bandwidth.  Another  common  complexity 
measure,  SampEn,  was  introduced  by 
Richman  and  Moorman  (33)  and  is  based 
on  ApEn  with  a  computational  refine¬ 
ment  (Appendix).  In  our  study,  ApEn  and 
SampEn  both  decreased  with  HS  and 
were  restored  with  resuscitation.  Previ¬ 
ously,  Palazzolo  et  al.  (28)  documented  a 
decrease  in  ApEn  secondary  to  hypoten¬ 
sion  in  dogs.  Using  selective  and  com¬ 
bined  blockade  of  the  autonomic  system 
branches,  the  authors  showed  that  under 
various  physiologic  conditions  (rest, 
standing,  hypotension),  entropy  reflects 
the  parasympathetic  modulation  of  the 
heart  rate  (28).  In  a  study  in  healthy 
humans,  Penttila  et  al.  (34)  found  that 
vagal  blockade  decreased  both  the  ApEn 
and  the  fractal  scaling  of  the  RRI.  Simi¬ 
larly,  in  our  study  the  decrease  in  various 
complexity  measures  was  associated  with 
evidence  of  vagal  withdrawal. 

The  fractal  dimension  of  the  RRI  time 
series  is  another  complexity  measure  ex¬ 
plored  in  this  study.  As  measured  by 
FDCL  (Appendix),  the  RRI  fractal  dimen¬ 
sion  changed  nonsignificantly,  although 
in  concert  with  the  other  complexity 
metrics.  These  results  are  similar  to  data 
presented  by  Kuusela  et  al.  (11),  in  which 
FDCL  decreased  along  with  ApEn  second¬ 
ary  to  terbutaline  administration  in  hu¬ 
mans.  Also,  West  et  al.  (35)  found  a  de¬ 
cline  in  RRI  fractal  dimension  in  subjects 
during  central  hypovolemia  induced  by 
LBNP. 

The  SymDyn  and  bits  per  word  meth¬ 
ods  used  here  for  quantifying  complexity 
are  computationally  quite  different  from 
the  ApEn  and  fractal  methods  (Appendix). 
Nevertheless,  changes  in  these  variables 
paralleled  those  of  the  entropy  and  fractal 
dimension.  Palazzolo  et  al.  (28)  showed 
decreased  SymDyn  during  rest,  standing, 
and  systemic  hypotension  in  dogs  as  a 
function  of  vagal  withdrawal.  SymDyn 
was  decreased  by  terbutaline  infusion  in 
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Figure  4.  A,  changes  in  R-to-R  interval  high-frequency  power  over  time.  Triangles,  hemorrhaged 
animals;  circles,  controls,  /-tests  (log-transformed  data):  *p  =  .004;  **p  <  .001;  ***/?  =  .015.  B, 
changes  in  systolic  arterial  pressure  (SAP)  high-frequency  power  over  time.  Triangles,  hemorrhaged 
animals;  circles,  controls,  /-tests  (log-transformed  data):  *p  =  .001;  **p  =  .019;  ***p  =  .012;  f p  = 
.004.  C,  changes  in  R-to-R  interval  integrated  baroreflex  sensitivity  over  time.  Triangles,  hemorrhaged 
animals;  circles,  controls.  Mann-Whitney  U  tests:  */?  =  .029;  **p  <  .001;  ***p  =  .002. 


humans  (11)  most  likely  due  to  decreases 
in  parasympathetic  cardiovascular  regu¬ 
lation  (36). 

In  contrast  to  the  RRI  metrics,  mea¬ 
sures  of  SAP  complexity  did  not  change 
in  our  study — yielding  results  similar  to 
those  of  Butler  et  al.  (3),  who  identified  a 
decrease  in  RRI  but  not  in  ABP  complex¬ 
ity  in  healthy  humans  during  LBNP  as¬ 
sessed  by  detrended  fluctuations.  As  sug¬ 
gested  by  Kuusela  et  al.  (11),  the  RRI  and 
SAP  may  have  different  underlying  dy¬ 
namic  structures,  with  the  former  being 
more  complex  and  the  latter  more  regu¬ 
lar  and  much  simpler  then  the  regulation 
of  the  heart  rate. 

Frequency-domain  HRV  analyses  have 
been  extensively  applied  to  critically  ill  hu¬ 
mans  and  have  yielded  important  findings. 
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Goldstein  et  al.  (7)  found  that  illness  sever¬ 
ity  scores  correlated  inversely  with  LFP  and 
HFP  of  the  RRI  in  critically  ill  and  injured 
pediatric  patients.  The  same  group  also  dif¬ 
ferentiated  brain-dead  and  non-brain-dead 
pediatric  patients  with  severe  brain  injury 
using  HRV  analysis  (37).  Winchell  and  Hoyt 
(8)  found  that  low  sympathetic  activity  rel¬ 
ative  to  vagal  activity  predicted  death  in 
ICU  patients,  whereas  sympathetic  predom¬ 
inance  was  associated  with  survival.  Cooke 
et  al.  (6)  applied  HRV  analysis  to  trauma 
patients  during  prehospital  transport. 
Death  was  predicted  by  parasympathetic 
predominance. 

In  contrast  with  the  previously  men¬ 
tioned  studies  in  which  relative  vagal  pre¬ 
dominance  predicted  mortality,  the  de¬ 
crease  in  RRI  HFP  observed  in  our 


nonlethal  study  suggests  that  progres¬ 
sive,  compensatory  cardiac  parasympa¬ 
thetic  withdrawal  occurred  in  response  to 
HS.  Previously,  Kawase  et  al.  (38)  found 
that  both  LFP  and  HFP  decreased  with 
progressive  blood  loss  in  anesthetized, 
mechanically  ventilated  dogs.  Similarly, 
Butler  et  al.  (3)  documented  cardiac  para¬ 
sympathetic  withdrawal  and  sympathetic 
stimulation  secondary  to  LBNP  in  hu¬ 
mans.  Triedman  et  al.  (39)  also  estab¬ 
lished  a  significant  decrease  in  vagal 
modulation  of  the  heart  in  humans  dur¬ 
ing  blood  donations. 

Lack  of  significance  in  the  changes  in 
LFP/HFP  ratio  in  our  study  may  reflect 
the  effects  of  anesthesia  on  HRV,  which 
alters  sympathetic  compensation  for  HS 
(5,  40,  41).  In  our  study,  integrated 
baroreflex  sensitivity  decreased  with  pro¬ 
gression  of  HS  and  was  reversed  with 
resuscitation.  The  changes  in  the  BRS 
that  we  observed  are  similar  to  those  re¬ 
ported  by  Iwasaki  et  al.  (42)  in  humans 
subjected  to  head-down  tilt,  in  which  a 
decrease  in  BRS  paralleled  the  decrease 
in  RRI  HFP.  Loss  of  HRV  complexity  in 
our  study  was  associated  with  a  decrease 
in  vagal  input  to  the  heart,  likely  a  nor¬ 
mal,  compensatory  response  to  hypovole¬ 
mia.  However,  we  cannot  exclude  the 
possibility  that  loss  of  complexity  at  lower 
blood  pressures  may  also  represent  isch¬ 
emia-induced  uncoupling  of  cardiovascu¬ 
lar  and  autonomic  regulation  (43)  and 
impairment  of  efferent  neuronal  input  to 
the  heart  and  lungs  (32,  44,  45).  Because 
brain  injury  has  also  been  characterized 
by  decomplexification  and  low  HF  and  LF 
powers  (7),  more  data  are  needed  to  en¬ 
able  discrimination  of  brain-injured  and 
HS  patients.  More  data  would  also  be 
needed  to  determine  the  sensitivity  and 
specificity  of  these  new  diagnostic  mea¬ 
sures. 

We  also  found  increases  in  SAP  HFP 
during  HS.  Julien  et  al.  (46)  suggested 
that  respiratory-frequency  oscillations 
(HFP)  in  arterial  pressure  involve  fluctu¬ 
ations  in  beat-to-beat  stroke  volume  of 
mechanical  origin.  ABP  oscillations  were 
found  to  be  mainly  caused  by  the  me¬ 
chanical  effect  of  respiration  on  intratho- 
racic  pressure.  As  shown  by  Saul  et  al. 
(47)  and  deBoer  et  al.  (48),  the  effect  of 
intrathoracic  pressure  on  stroke  volume, 
in  turn,  is  enhanced  by  hypovolemia.  Au¬ 
tonomic  regulation  of  the  HR  may  also 
cause  vagally  mediated  changes  in  the 
SAP  HFP  (47). 

As  shown  by  Brown  et  al.  (29),  breath¬ 
ing  frequency  and  depth  are  known  to 
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influence  HRV  profoundly  and  should  be 
controlled  in  HRV  studies.  Because  our 
animals  were  breathing  at  a  set  rate  and 
depth,  the  influence  of  lung  inflation  on 
RRI  responses  was  constant  throughout 
the  experiment.  Thus,  we  can  attribute 
the  changes  in  complexity  in  this  exper¬ 
iment  to  autonomic  rather  than  mechan¬ 
ical  sources  (30).  Changes  in  heart  rate 
complexity  might  be  different  in  patients 
who  are  not  mechanically  ventilated. 

A  limiting  factor  affecting  HRV  met¬ 
rics  is  the  nonstationarity  of  the  signal. 
For  a  stationary  time  series,  the  average 
value  of  the  mean  and  the  sd  do  not 
vary  with  time  (15).  Signals  that  do  not 
obey  this  condition  are  considered  non¬ 
stationary. 

Nonstationarity  affects  frequency-do- 
main  variables  more  than  it  does  nonlinear 
variables.  The  latter  are  not  immune  to  the 
effects  of  signal  nonstationarity  (especially 
fractal  scaling).  However,  some  of  them — 
especially  the  entropy-based  measures — 
quantify  statistical  order,  as  they  were 
designed  for  quantification  of  regularity 
in  both  deterministic  and  stochastic 
models  (15).  New  methods  not  explored 
in  the  present  study,  such  as  detrended 
fluctuation  analysis  that  quantifies 
changes  in  short-  and  long-scale  compo¬ 
nents  of  data  sets  (49),  may  be  more  ad¬ 
vantageous  for  analysis  of  nonstationary 
signals. 

In  an  effort  to  minimize  the  influence 
of  nonstationarity  on  the  data  sets  ana¬ 
lyzed  in  this  study,  instead  of  continuous 
data  sampling  during  protocol  interven¬ 
tions  we  chose  to  analyze  waveforms  re¬ 
trieved  from  periods  of  time  between  in¬ 
terventions  during  which  the  data  were 
as  stable  as  possible.  However,  the  poten¬ 
tial  effect  of  signal  nonstationarity  on  our 
results  cannot  be  discounted,  and  our 
data  should  be  interpreted  with  caution. 

The  presence  of  ectopic  beats  may  al¬ 
ter  the  results  of  HRV  analysis.  Our  high¬ 
est  rates  of  interpolation  were  1.5%  and 
1.8%  (in  two  data  sets  of  120  analyzed); 
thus,  we  believe  that  the  minimal  edit¬ 
ing  of  the  ectopic  beats  performed  in 
our  study  did  not  significantly  affect 
our  results. 

Currently,  techniques  are  sought  that 
enable  early  identification  of  subjects 
with  “survivable”  or  “nonsurvivable”  in¬ 
juries,  enabling  distant  decision  support 
and  increased  situational  awareness  in 
mass-casualty  or  battlefield  settings,  a 
concept  known  as  remote  triage.  Assess¬ 
ment  of  HRV  with  nonlinear  statistical 
methods  may  be  one  such  technique.  By 
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providing  insight  into  the  underlying 
control  mechanisms  at  work  during  hem¬ 
orrhagic  shock,  it  is  possible  that  HRV 
analysis  may  augment  the  information 
obtained  from  traditional  vital  signs  such 
as  the  pulse  and  blood  pressure.  However, 
studies  in  human  trauma  patients  are 
needed  and  are  ongoing.  Continuous  data 
analysis  with  sliding  time  windows  for 
both  spectral  and  nonlinear  analysis  (not 
employed  in  our  study)  may  help  clarify 
whether  and  at  what  point  HRV  variables 
may  be  superior  to  traditional  vital  signs. 

CONCLUSIONS 

This  study  in  anesthetized  swine  sup¬ 
ports  the  concept  that  RRI  complexity  is 
decreased  by  hemorrhagic  shock.  Several 
complementary  measures  of  complexity 
were  noted  to  change  with  shock,  to  in¬ 
clude  ApEn,  SampEn,  and  SymDyn. 
Changes  in  FDCL  were  similar  but  not 
significant.  These  alterations  in  the  com¬ 
plexity  of  the  RRI  may  be  mediated  by 
vagal  withdrawal,  as  evidenced  by  con¬ 
comitant  decreases  in  RRI  HFP  and  in 
integrated  baroreflex  sensitivity.  These 
techniques  are  now  being  applied  to  the 
analysis  of  heart  rate  variability  in  trauma 
and  bum  patients. 
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APPENDIX 

Baroreflex  Sensitivity  Calculation 

Baroreflex  sensitivity  is  calculated  au¬ 
tomatically  by  identification  of  a  mono- 
tonically  increasing  or  decreasing  set  of 
three  or  more  consecutive  corresponding 
data  points  on  the  RRI  and  SAP  time 
series.  From  the  set  of  SAP  and  RRI  val¬ 
ues,  an  xy-scatter  plot  is  created  and  a 
regression  line  is  fitted  over  the  data 
points.  The  slope  of  this  line  is  the  BRS. 

Nonlinear  Variables  (11) 

Approximate  Entropy.  ApEn  is  a  mea¬ 
sure  of  the  complexity  of  a  signal;  it  is  low 
in  regular  signals  and  high  in  irregular 
signals.  ApEn  determines  the  conditional 


probability  of  finding  specific  patterns  in 
the  time  series,  i.e.,  the  logarithmic  likeli¬ 
hood  that  a  run  of  patterns  that  is  close 
remains  close  on  the  next  incremental 
comparison.  The  template  patterns  are 
constructed  from  the  signal  itself,  and  no  a 
priori  knowledge  of  the  system  is  needed. 

Sample  Entropy.  This  is  a  similar  con¬ 
cept  to  ApEn,  with  the  computational  dif¬ 
ference  that  the  vector  comparison  with 
itself  is  removed.  ApEn  occasionally 
falsely  indicates  that  signals  are  more 
regular  than  they  actually  are.  SampEn 
corrects  for  this  problem. 

Symbol  Dynamics  Entropy  (Sym- 
Dyn).  Four  “symbols”  are  generated  (A,  B, 
C,  D)  by  using  the  mean  and  sd  of  the 
signal  to  create  value  “bands”  for  each  sym¬ 
bol:  e.g.,  “A”  if  signal  value  <  (mean  -  SD), 
etc.  The  time  series  is  then  converted  to  a 
symbol  sequence.  This  symbol  sequence  is 
then  interpreted  as  a  sequence  of  three- 
letter  “words.”  The  complexity  (Shannon 
entropy)  of  the  distribution  of  words  is  then 
determined. 

Bits  per  Word.  Bit  per  word  is  a 
method  of  normalizing  SymDyn. 

Lempel-Ziv  Entropy.  The  signal  is 
converted  to  a  binary  sequence  (series  of 
0s  and  Is),  in  which  1  represents  values 
above  the  mean  and  0  represents  values 
below  the  mean.  This  kind  of  sequence 
can  be  constructed  by  two  operations:  the 
insertion  of  symbol(s)  to  form  a  subqueue 
(i.e.,  a  shorter  sequence  of  0s  and  Is)  and 
the  copying,  or  repetition,  of  this  sub¬ 
queue.  For  Lempel-Ziv  entropy,  the  com¬ 
plexity  of  this  sequence  is  determined  by 
the  number  of  such  insertion  and  copy¬ 
ing  operations  that  need  to  be  performed 
to  generate  the  sequence. 

Fractal  Dimension  by  Dispersion 
Analysis  (FDDA).  FDDA  determines  the 
fractal  dimension  of  a  time  series.  A  new 
time  series  is  created,  based  on  the  mean 
of  each  set  of  two  adjacent  values.  This 
process  is  continued  with  data  sets  of  four 
points,  16  points,  etc.  Log  SD(m)  is  plot¬ 
ted  against  log  m.  FDDA  is  defined  as 
1-slope.  The  variable  m  is  the  number  of 
data  points  used  to  group  adjacent  data 
points:  4,  16,  etc.  FDDA  has  a  value  of 
1.0— 1.5.  FDDA  =  1  indicates  a  constant 
signal,  and  1.5  indicates  a  maximally  frac¬ 
tal  structure,  i.e.,  a  totally  random  signal. 

Fractal  Dimension  by  Curve  Lengths. 
Another  method  of  determining  the  frac¬ 
tal  dimension  of  a  time  series.  A  given 
curve  is  conceptualized  as  consisting  of  a 
number  of  short  segments.  FDCL  counts 
the  number  of  such  segments  of  various 
lengths  needed  to  follow  the  curve. 
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